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We report a theoretical study of time-dependent transport in a ballistic graphene field effect 
transistor. We develop a model based on Floquet theory describing Dirac electron transmission 
through a harmonically driven potential barrier. Photon-assisted tunneling results in excitation of 
quasibound states at the barrier. Under resonance condition, the excitation of the quasibound states 
leads to promotion of higher-order sidebands and enhanced higher harmonics of the source-drain 
conductance. The resonances in the main transmission channel are of the Fano form, while they 
are of the Breit-Wigner form for sidebands. We discuss the possibility of utilizing the resonances in 
prospective ballistic high-frequency devices, in particular frequency multipliers. 


PACS numbers: 72.80.Vp, 73.23.-b, 73.23.Ad 

Already in the early years of graphene research, ana¬ 
logue high-frequency electronics was recognized as a po¬ 
tential niche for applications m- Although many 
devices have probably been limited by parasitics due 
to problems with developing good recipes for making 
graphene transistors, the current speed record [5] is al¬ 
ready a cut-off frequency of over 400 GHz. At the same 
time, we have seen a rapid improvement of graphene ma¬ 
terial quality. Mobilities reaching 10 5 cm 2 /V-s at room 
temperature and larger than 10 6 cm 2 /V-s at low temper¬ 
ature have been achieved |4[. Promising paths towards 
improved mobility include encapsulation of graphene be¬ 
tween layers of other two-dimensional (2D) crystals, no¬ 
tably hexagonal boron-nitride, or suspension of graphene 
between contacts. Using the latter approach, ultra high- 
quality p — n junctions were recently made [6|. Fabry- 
Perot resonances at zero magnetic field were measured, 
and so-called snake states were possible to see at small 
magnetic fields of order 20 mT. With such rapid improve¬ 
ments of device quality, it has become increasingly impor¬ 
tant to study in detail ballistic high-frequency devices. 

One of the key ideas behind using 2D materials for 
high-frequency electronics is the favorable scaling to¬ 
wards short gate lengths without so-called short-channel 
effects pQ. Thin channels (2D is the extreme) allows 
for short gates, high speed, and high-density integration. 
High speed, reaching THz frequencies [ 7 ], is the ultimate 
goal. Another advantage of graphene is the possibility 
to tune the electron density, for instance by means of a 
back gate: the Fermi energy can be tuned from the elec¬ 
tron to the hole band (through the so-called Dirac point 
at charge neutrality). Such ambipolarity is very advan¬ 
tageous, in that both n-type and p-type devices can in 
principle be made at will across a single wafer. 

A challenge is to capitalize on the unique properties of 
graphene and derive device functionality directly from 
the fact that electrons in graphene behave like mass¬ 
less Dirac particles with linear spectrum and a pseu¬ 
dospin degree of freedom. Several works in this direc¬ 
tion show indeed that ac transport in graphene is a 


rich subject. Studies include quantum pumping isma, 
non-linear electromagnetic response [13HT8] , and photon- 
assisted tunneling phenomena jl9ti25] . In theoretical in¬ 
vestigations for low doping (Fermi energy Ep close to 
the Dirac point) and high frequencies £2, with Ep and Cl 
of comparable magnitude (we put h = 1), a true quan¬ 
tum mechanical description becomes necessary. For time- 
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FIG. 1. (a) A graphene field effect transistor, where the 

overall doping level is controlled by a back gate (BG), and 
the source (S) - drain (D) current is controlled by the top 
gate (TG) dc and ac signals, (b) The harmonic ac signal 
of frequency U leads to inelastic scattering that under reso¬ 
nance condition excites an otherwise unoccupied bound state 
in the top gate barrier potential at energy Eb- This leads to a 
Fano resonance in transmission to Eq due to interference be¬ 
tween processes (1) and (2) and a Breit-Wigner resonance in 
transmission to E 2 [process (4)]. Process (4) leads to higher- 
harmonic generation, in particular the 2Q harmonic. 
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dependent transport in two-dimensional electron gases in 
semiconducting heterostructures, displaying a quadratic 
dispersion relation, photon-assisted tunneling in time- 
harmonic potentials is described well within a Floquet 
theory framework and have been investigated for a long 
time [26l[28l . Here, we study theoretically a ballistic field 
effect transistor with a harmonic drive applied to the top 
gate, see Fig. [lja), within a Floquet theory applicable 
to graphene. The harmonic drive of frequency D sup¬ 
ports inelastic scattering from the Fermi energy Ep, to 
sideband energies E n = Ep + ttD, where n is an inte¬ 
ger. Near charge neutrality, on the scale of the drive fre¬ 
quency, the barrier is close to transparent due to Klein 
tunneling. At the same time, a quasibound state on the 
barrier can be inelastically excited through a resonant 
process (supported by the harmonic drive) that interfere 
with direct elastic transmission. This leads to a Fano res¬ 
onance in direct transmission, as recently found numeri¬ 
cally [23M25] . Here, we show that higher order sidebands 
are simultaneously resonantly enhanced, which leads to 
the possibility of building a frequency multiplier based 
on a ballistic graphene device that we study in detail in 
this paper. 

We are interested in the intrinsic properties of the 
graphene transistor in Fig. [lja), and neglect parasitics. 
This allows us to make a minimal model in terms of a 
Dirac Hamiltonian 

H = ~icr x V x + cTyky + [Z 0 + Z\ cos(f2t)] S(x), (1) 


the Landauer-Biittiker approach [29] is used to compute 
the current operator in terms of creation and annihila¬ 
tion operators for incoming and scattered waves, where 
the latter are related to the former through the scattering 
matrix. A statistical average is performed to obtain the 
time-dependent current that depends on the occupation 
factors of the source and drain leads, which are given by 
the Fermi function with chemical potentials shifted by 
the applied voltage eV. Complete derivations of all for¬ 
mulas are given in the Supplemental Material [30] . Be¬ 
low, we shall give results for conductance in linear re¬ 
sponse to the applied voltage at zero temperature. 

For the scattering problem, since the Hamilto¬ 
nian is periodic in time, we use a general Floquet 
ansatz i/;(x,k y ,i) = 'ip n (x, k y , E) exp(— iE n t). When 
plugged into the Dirac equation it yields a set of differen¬ 
tial equations for the (formally infinitely many) sideband 
amplitudes k y ,E). In the following we do not write 

the arguments x, k y , and E in order to keep the notation 
compact. The sideband amplitudes can be arranged into 
a vector <f> = [..., Vh? ...] T , which then satisfies 

= M t d<&, where 

Z 

M-td — \k , y&z T ^D n G x ^ZoG x 8{Ej\ (S) 1 i 5{E)(j x (S) 2 

is a tridiagonal matrix in sideband space with (l) nm = 
3nm £md (2) nm — ^n,m-\-i T ^n,m —i • After integration over 
x = 0 we obtain a boundary condition 


where we have set the Fermi velocity in graphene equal to 
unity, vp = 1. The Pauli matrices are as usual denoted 
<j x and a y . The top gate barrier potential is considered 
smooth on the atomic scale and cannot induce scattering 
between the two valleys in the band structure. In the 
end all observables will contain an extra factor of two 
to account for valley degeneracy, in addition to spin de¬ 
generacy. At the same time, on the Dirac length scale 
(given by hvp/Ep after reinstating h and vp), we con¬ 
sider the potential width D to be small but its height 
V to be large, such that we can take the limits D 0 
and V oo keeping the product VD = Z constant. The 
strengths of the time-independent component Zo and the 
time-dependent component Zi can be different. The 5- 
function in Eq. 0 is therefore smooth on the atomic 
scale but sharp on the Dirac length scale. We consider 
the barrier to be translational invariant along the trans¬ 
verse direction, which guarantees that the corresponding 
wave vector component k y is conserved. The spatial de¬ 
pendence then enters through the coordinate x perpen¬ 
dicular to the barrier. Finally, we assume homogeneous 
doping of the graphene sheet tuned by the back gate. 

The methodology to solve the problem at hand is 
to first solve the scattering problem for the wave- 
functions satisfying the Dirac equation l~Lil;(x,k y ,t) = 
id t 'ip(x,k y ,t). The solution can be collected into a uni¬ 
tary scattering matrix for reflection and transmission co¬ 
efficients between incoming waves at energy E and scat¬ 
tered waves at energies E n = E + nVt. For the current, 


&(x = 0 ) = exp 


iZoG x (£) 1 + i 



$(x = 0 + ) 


= MT(x = 0 + ). 


( 2 ) 


This boundary condition can also be derived by solving 
a square barrier problem first, and in the end let D —>> 0 
and V oo keeping the product VD = Z constant, see 
also Ref. [31]. This boundary condition gives an elegant 
view of scattering off a potential in graphene in terms of 
pseudospin rotation. For instance, the static barrier leads 
to a rotation around the pseudospin x-axis by an angle 
— 2Zq. If the pseudospin is aligned with a x , i.e. elec¬ 
tron propagation along the x-axis with perpendicular in¬ 
cidence, the rotation has no effect (Klein tunneling [32]). 
For other angles, transmission is non-perfect. The ex¬ 
plicit formulas for the transmission amplitudes t n (k y ,E) 
derived from Eq. Q are given in Ref. [30] . 

In Fig. [2] we display the transmission probabilities 
T n (E,(p) = \t n (E, (p)\ 2 for n = 0 and n = 2, where 
T n (E,ip) denotes incidence on the barrier at energy E 
and transmission at sideband energy E n , keeping the 
parallel momentum k y = |£7| sin conserved (the an¬ 
gle (f is measured relative to the barrier normal). In 
the main transmission channel To(E,ip), Klein tunnel¬ 
ing is apparent in that the transmission is very close to 
unity. Deviation from unity transmission is due to the 
static barrier of strength Zo and finite incidence angle 
(non-zero ip) and, in addition, scattering to other side¬ 
bands with n / 0. The transmission probability to the 
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E r = E b -f f2, see Fig. |2jb). The resonance in T 2 (E,p) 
can be viewed as due to transmission in energy space 
through a double barrier structure with barrier heights 
proportional to Z\. 

To extract more information about the above numer¬ 
ical results, we proceed with an analytic analysis. We 
can expand the boundary condition in Eq. © to second 
order in the ac drive strength Zi, assuming Zi <C 1, 


M 


p iZQ<j x 


T _L *'^ 1 
1 + ~2 ax 


z 2 

—p(2 • I + 3) 


(3) 


where (3) nm = £ n?m+ 2 + 5 n ,m- 2 in sideband space. To 
second order in Z i, the transmissions to the first two 
sidebands can be computed Go] by solving a system of 
equations for to, t± i, and t± 2 . We separate two cases: (i) 
off-resonant transmission and (ii) on-resonant transmis¬ 
sion. For case (i) off-resonant transmission, the equation 
system can be inverted directly and we get (for each p; 
we suppress the argument p below for brevity), 


to(E) 


1 + ^ + Z^°\E)Ao, 1 (E)t^(E 1 )A h0 (E) 
+t<- 0 \E)A o ,- 1 (E)tW(E- 1 )A- lt0 (E)] tW(E), 


FIG. 2. Energy and incidence angle dependence of transmis¬ 
sion probabilities for (a) elastic scattering To(E, ip) and (b) in¬ 
elastic scattering between energy E and E+ 2D, T 2 (E, p). The 
black region to the left of the white dashed lines in (b) are re¬ 
gions where the second sideband wavefunctions are evanescent 
waves decaying away from the barrier. The barrier strengths 
are Zq — 0.47T and Z\ — 0.45. Inset: transmission probabili¬ 
ties for fixed p — 7r/9. 


second sideband T 2 (E,p) is in general very small. For 
certain energies there are Fano resonances |23ff25l in¬ 
duced by the time-dependent drive and a bound state 
at the barrier, which give rise to a peak-dip structure 
dispersing with p, one feature at positive energies and 
another one at negative energies. The Fano resonances 
occur in a parameter range where the outgoing (from 
the barrier) wavefunctions at a sideband E n are evanes¬ 
cent (below we shall concentrate on n = ±1, which are 
the most pronounced resonances in Fig. [2]). This hap¬ 
pens when inelastic scattering from E = k x (E ) 2 + k 2 to 
E n = k x (E n ) 2 + k 2 (with conserved k y ) causes k x (E n ) 
to become imaginary. Resonant behavior occurs due to 
the existence of a bound state on the barrier at energy 
Eb(Zo,ky) = — sgn(Zo)\ky\ cosZ q, see Ref. [301, that can 
be excited by the ac drive (in which case it becomes quasi¬ 
bound). The Fano resonance at E r = E^Efl is a quantum 
mechanical interference between direct elastic tunneling 
and a tunneling process involving excitation to the first 
sideband (for n = =L1 at E r = E^ =p D) and de-excitation 
back to energy E, see paths (1) and (2) in the diagram 
in Fig. [ljb) for the E < 0 case. On resonance, inelastic 
tunneling to the second sideband is resonantly enhanced 
and T± 2 (E, p) display a Breit-Wigner resonance peak at 


t ±l {E) « -ZP°\E ±1 )A ±lfi (E)t^(E), 

t ±2 (E) » Z\ [t^(E ±2 )A ±2 , ±1 (E)t^(E ±1 )A ±lfi (E) 

-t^(E ±2 )A ±2fi (E)]t^\E), (4) 


where t^\E n ) is the transmission amplitude without ac 
drive computed at energy E ni and A njm (E) is a transi¬ 
tion amplitude in energy space between energies and 
E ni which can be related to off-diagonal matrix elements 
(in sideband space) of the matrix M in Eq. ©. The 
above expressions make the inelastic tunneling processes 
at play explicit, see enumerated processes in Fig. 0b). 
For instance, the expression for t\{E) read from right 
to left has a transparent physical meaning. It consists 
of transmission amplitudes at E and E\, separated by 
a transition in energy space Ai : 0 , corresponding to ab¬ 
sorption of one quantum Q. Consequently, the process is 
of order Z\. Direct transmission has corrections to the 
static transmission amplitude due to excitation and de¬ 
excitation to neighboring sidebands (processes of order 
Z\\ while t 2 (E) consists of a direct process of absorbing 
two quanta, 2D, and a sequential process involving the 
first sideband energy, both are of order Z\. This tells us 
that the sideband amplitudes are in general very small 
when Z\ is small. 

The above picture changes for case (ii) on-resonant 
transmission, for energies near E r = E^ =L D [we shall 
concentrate on E^ — D in the following discussion, as in 
Fig. [ljb)]. In this case, the equation determining the 
function t^(E n ,k y ) at energy E n = E & (here n = 1), 
has to be reconsidered. There is a pole in the matrix 
equation determining the scattering matrix at this en¬ 
ergy for fixed k y , corresponding to formation of a bound 
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state with evanescent waves decaying away from the bar¬ 
rier. The bound state is unoccupied (decoupled from 
reservoirs) in the absence of ac drive. For case (ii) on- 
resonant transmission, we get for energies SE around the 
resonance energy E r 


to(E r + 5E ) 
£2 {E r + 5E ) 


^ ~ Zl^jEr) (Q) , p . 

5E + Z*hi(E r ) 1 b> ’ 


ZjhsjEr) 

SE + Zfhi(E r ) 


t(°\E b + 2n), 


(5) 


where hi(E r ) = hi(E r ,k y , Zq), i = 1,2,3, are complex 
functions given in Ref. [2Q] (their explicit form is not 
important in the discussion below). Note that t\ is not 
well defined near resonance (it was eliminated in the cal¬ 
culation) because it is related to the excitation of the 
bound state. The conductance computed below will not 
get contributions from this sideband energy [crossed pro¬ 
cess (3) in Fig. [ljb)]. For the direct transmission prob¬ 
ability To(E r + SE ), neglecting for a while the second 
sideband contribution (setting /i 2 = 0 above), there is 
a characteristic Fano resonance form To(E r + SE) oc 
(qT / 2 + SE)/[SE 2 + (T/2) 2 ], where T oc Z 2 and q is of 
order unity, O [(Z 1 ) 0 ]. This is the blue dotted line dis¬ 
played in the inset of Fig.[2][b). Taking into account tun¬ 
neling (in energy space) to the second sideband (h 2 finite 
above) and higher, we obtain the corrected line-shape, 
the black solid line in the inset of Fig. [2](b) . For the prob¬ 
ability to scatter inelastically to the second sideband, we 
obtain from above a Breit-Wigner resonance with the 
characteristic form T 2 (E r + SE) oc (T/2) /[SE 2 + (T/2) 2 ], 
which is displayed as the black dashed line in the inset 
of Fig.^b). Thus, in a range SE oc Z 2 Tl around E r , the 
response is highly non-linear and higher-order harmonics 
can be resonantly enhanced. 

To quantify this, we present calculations of the lin¬ 
ear conductances G n , both the time-averaged compo¬ 
nent (n = 0) and the first harmonics (n > 0), see 
Fig. HI Note that in linear response (small source-drain 
voltage), the source drain ac current / = ^ n / n e _mm , 
with its harmonics / n , naturally define ac conductance 
components G n , see formulas and additional figures in 
Ref. [30]. In Fig. [3](a) we plot the angle resolved conduc¬ 
tance Go(Ef,(p), which reflects the sum over transmis¬ 
sion functions in Fig. |2j After angle integration, the dc 
linear conductance (thin green straight lines) changes due 
to the ac drive into the solid black line in Fig. |3](b) . The 
Fano resonance is clearly visible as a peak-dip feature in 
Go(Ep). Thus, it is enough to study the time-averaged 
conductance to infer influence of the ac drive. In Fig. [3](b) 
we also present the real and imaginary parts of the sec¬ 
ond harmonic G 2 (£ , j p). For small drive amplitude Zi, 
the harmonics generally scale as Zf in perturbation the¬ 
ory and the second harmonic is expected to be small. 
Near resonance, however, it is enhanced to order unity, 
O [(Z 1 ) 0 ], within a window of doping ^ Z 2 Q around 
Ep = E r . The results in Fig. [3] are obtained for a rather 


stronger drive Zt =0.45, including many sidebands. 
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FIG. 3. Source-drain linear conductances in the presence of ac 
drive on the top gate. Upper panel: impact angle resolved av¬ 
erage conductance Go(Ef,<p)- Lower panel: angle integrated 
average conductance Go(Ef) and real and imaginary parts of 
the second harmonic G 2 (Ff) (multiplied by a factor 10). The 
dip-peak structures are related to the Fano and Breit-Wigner 
resonances in the elastic and inelastic transmission functions. 
Near resonance, the second harmonic is enhanced from 0[Zf] 
to order unity, where Zi is the ac drive strength. 


the enhancement of G 2 is visible in a wide range of ener¬ 
gies. We note that there is an unfavorable prefactor for 
the ac components, reflected in the transmission function 
T 2 (E) being of order 0.2 even on resonance, see Fig. [2](b) 
and formulas in Ref. [3Qj. 

In summary, we have investigated time-dependent 
transport in a ballistic graphene field effect transistor 
with ac drive on its top gate. We find resonances in 
inelastic scattering to sideband energies, related to exci¬ 
tation of a quasi-bound state in the top gate barrier. This 
leads to substantial resonant enhancement of higher har¬ 
monics of the source-drain conductance, that could pos¬ 
sibly be used in developing a frequency multiplier based 
on a ballistic device. 
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Supplemental Material 

I. WAVE SOLUTIONS IN GRAPHENE 
A. General solution 

We start by introducing general wave solutions in graphene without time-dependent perturbation. They are known 
(see e.g. Refs. [#© and we write them down here to establish a coherent notation for subsequent sections. As 
mentioned in the main text, we consider only one valley (one A-point) described by the Hamiltonian 

Ho = -icr ■ V, a = (a x ,cry). (SI) 

We have to solve the Dirac equation idttp(x,y,t) = Hotp(x, y, t), which is done by the standard ansatz 

ip(x, y, t ) oc e ikxX e ikyV e~ iEt ip(k x , k y ,E). (S2) 

We obtain the following eigenvalues and eigensolutions, 

E x (k x ,k v ) = Xy/kl + k^, A = ±1, (S3) 

Mk x ,k y ,E) = —j= ^ kx+iky ^ . (S4) 


B. Scattering basis 


Note that once we have found the spectrum, Eq. (S3), there are only two independent parameters labeling eigen¬ 
states, e.g. (k x ,k y ) or (k y ,E). Since we are going to build a scattering theory following Biittiker [^3]-^5] the latter 
choice is natural because we assume translational invariance along the barrier (y-axis), c.f Eq. (1) in the main text. 
In order to introduce the scattering basis we have to find the group velocity of states propagating along the x-axis 
(perpendicular to the barrier). Using standard definitions we have, 


v(ky,E) = = ±v(ky,E), y(ky,E) = Kx ( k y’ E ) ' 


E 


(k v ,E)= S gn(E)^E^ki, 


(S5) 


where the upper and lower signs describe particles moving in the positive and negative directions along x, respectively. 
Then we can introduce a scattering basis via 


ky,E) = 


ip^(x,k y ,E) = 


1 


y/2v(ky,E) \v(ky,E) 


1 


y/2v(ky,E) \jl(ky,E) 


v(k y ,E) = 


K'xi.kyi E ) “l - iky 

E 


, r){ky,E) = 


in x (ky,E)x 

Ex (J^y , + iky 

E 5 


(S6) 


where arrows indicate the direction of propagation. The normalization in Eq.(S6) is chosen such that a particle 
incident from the left carries unit flux, defined as 


j x (x,ky,E) = ift(x,ky,E)a x ip(x,ky,E). 

This basis is used to find a scattering matrix and build the scattering field theory below. 


(S7) 


II. FLOQUET SCATTERING MATRIX IN GRAPHENE WITH AC 5-POTENTIAL 

Let us now discuss the Floquet scattering matrix for graphene in presence of an oscillating line scatterer, i.e. we 
consider a system described by [c.f. Eq. (1) in the main text] 

n = n 0 + [Z 0 + Z 1 cos (Sit)] 5(x). (S8) 

Before discussing the solution associated to the full time-dependent Hamiltonian, it is instructive to consider Zi = 0. 
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A. Static ^-barrier 

In this case, since scattering is elastic, it is easy to write down a scattering ansatz, assuming incoming particles 
from the left, 

J,( T t F] _ f^(x,ky,E)+A°'l^(x,k y ,E) ifx<0, . . 

n ' y ’ ’ ~ k y , E) if X > 0. (b9) 

The superscript indicates functions X computed for a static barrier. The unknown transmission t^°\k y ,E) and 
reflection A°\k y ,E) coefficients are found through the boundary condition at x = 0, which reads [c.f. Eq.(3) in the 
main text] 


^(0 , ky , E) = exp[iZ o cr ;E ]'0(O + , k y , E ). 


(S10) 


It is straightforward to find a solution to Eq.( |S10| ), but it is convenient for what follows to write down an equation 
satisfied by t^°\ 


D(ky,E)t^ = 1, 


D{ky ’ E) = M^E) !) eMiZo* x ] 


Using Eq.(S6) we can easily simplify Eq.(|Sll[) and obtain 


t(k y , E) = D(k y , E) 1 = ( cos Zq + i 


sinZ 0 

J v(k y ,E) 


-1 


If we introduce an incidence angle <p via k y = \E\ simp, then Eq.(S12) can be rewritten as 


f (0) (ky, E) = 


COS (p 


cos p cos Zo + i sin Zo 


(Sll) 


(S12) 


(S13) 


B. Barrier-induced bound state 


It is well-known [^6] that poles of the scattering matrix correspond to bound states. In our case the s tatic (5-barrier 
induces exactly one bound state as will be shown now. We equate to zero the denominator of Eq.(S12) and impose a 
condition that the bound state solution has to be decaying away from the barrier, which means 


(ky j T/) — iyjky E ^, Kj X (k y , E 1 ) — iE tan Zq . 


(S14) 


One can see that Eq.(S14) is periodic in Zo and we consider for definiteness —§ < Zo < f. Then the energy of the 
bound state is given by 


E b = —sgn(Z 0 )|fcy| cos Z 0 . (S15) 

It is interesting to note that the bound state plays no role in DC transport since it is disconnected from the continuum 
of propagating waves. This circumstance changes as soon as we allow inelastic scattering on the barrier, when Zi / 0. 


C. Oscillating (5-barrier 


In the case when Z i ^0 the Hamiltonian, Eq. (S8), is periodic in time, which enables us to use the Floquet theorem 
EJd7ld8] for finding eigenvectors, 


+oo 

ip(x, ky,t) = e~ lEt y] e~ tnSlt ip n (x,ky,E). 


n= — oo 


(S16) 













Now if we introduce a column vector 


$(x, ky, E) = (... ky, E),ipo(x, ky, E),ipi(x, k y , E ),... ) T , 


(S17) 


then the condition to be satisfied at x = 0 is [see also Eq.(3) of the main text], 

$(0-,ky,E)^M<f>(0 + ,ky,E), 


M = exp 


Z x 

iZ 0 cr x <g> 1 + i—(7x 


(S18) 


P] n ,ra _ P] n ,m 3n,m +1 +^n,m-l- 


Since the barrier is active only at x = 0 , asymptotic solutions are still given by a linear combination of the static solu¬ 
tions, Eq.(|S 6 |). The barrier only scatters an incident particle with quantum numbers (E, k y ) into a linear combination 
of states with quantum numbers ( E n ,k y ), where E n = E + (in the end we have to consider only propagating 
outgoing waves, E n > | k y \, for calculating transport properties). Therefore we use the following ansatz 


7p n (x,k y ,E) = 


S ni o^(x,k y ,E n ) +r n 'ip^(x,k y ,E n ) if x < 0, 

^ 71 * 0 —>■(*■£ 5 ky , E n ^j if X > 0. 


(S19) 


We can eliminate reflection coefficients r n (k y ,E ) and find a system of equations for t n (k y ,E ) only, which reads 


^ 2y/y(ky,E n )y(ky,E m ) ^ v{ky,En) X ) [ M L™ [ v (k y , E n 


trn ^n, 0 * 


(S 20 ) 


Eq. (S20) must be solved numerically. 


III. ANALYSIS OF SIDE-BAND TRANSMISSION COEFFICIENTS: FANO AND BREIT-WIGNER 

RESONANCES 


The system of equations ( |S20| ) is in principle infinite in sideband index space. To find an approximate solution we 
have to cut the system by setting a maximum allowed n max number of side-band. We assume that Zi <C 1, expand 
M up to terms of order O(Zf), and consider five outgoing channels with n = {0, ±1, ± 2 }. Then we obtain a system 
of five coupled equations which reads [omitting the arguments ( k y ,E ) for brevity] 

(l — D 2 t2 + Z 1 A 2 ,\t\ + Z\A 2 $t§ = 0, 

Z\A\^ 2 t 2 + ^1 — Diti + ZiAi ?0 to + Z^Ai-it-i = 0, 

< ZiAo, 2 t 2 + ZiAq^\ t\ + ^1 —Doto + ZiAo_it-i + ^ 1 ^- 0 ,- 2^-2 = 1 5 (S21) 

Z\A^i t iti + ZiA_i ? oto + ^1 — 4 ^ D-it-i + Z\A—\ — 2 t— 2 = 0, 

Z\A- 2 $t§ + ^1^-2,-i^-i + (l — D- 2 t -2 = 0, 

where we have used the following notations 


D n (ky, E) = D(k y , E n ), 


d — 


(i/2)l n m l 


n — ml 


, 7 7 T F \ ft t-i T (~v(ky’ E n) 1 ) exptiZoggpl." m| ( / fc ^ b, n ^ m. (S22) 
j y^yky , E n j v(k y , E m ) y> m) J 


Note that from Eq. (S 22 ) and Eq. (Sll) it is obvious that D“ 1 (/c ?/ ,E) = t(°\k y ,E n ) provided the corresponding 
wave is propagating, i.e. E n > \k y \. On the other hand the new functions A n , m have a meaning of transition matrix 
between the side-bands. Now we recall that the presence of a (static) ^-barrier implies existence of a bound state, 
see Sec |IIB[ which now can be coupled to the propagating waves via inelastic scattering. In this case one of the 
functions D n (k y , E) vanishes when E n = E^. This possibility leads to resonances in the transmission spectrum of the 
side-bands [see Fig|Sl|, as will be discussed in details below. 
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FIG. SI. Energy and incidence angle dependence 


of transmission probabilities for inelastic scattering to the side-bands. 


A. Off-resonant transmission 


We will first consider the rather trivial case of transmission in different side-bands away from the resonances. In 
this case we can straightforwardly estimate orders of magnitude for the side-band transmission coefficients keeping 
only contributions 0(Z\\ 


to = 4 0) + Z 1 T 0, 
t±l = Z\T±\, 
t ± 2 = Zfr± 2 , 


(S23) 


where we introduced for convenience tn\k y ,E) = t^\k y ,E n ). Keeping the same order of approximation in Eq.(S21) 
we can easily solve it with the following result 


r 0 =(I+C^o,i 


r±i = - 

r± 2 = 


(°) 4- . + +(°) G- 1 ’ 0 ^ t (0) 

D 1 +t ° °’- 1 ) ° ’ 

^±l, 0 + ( 0 ) 

( 0 ) 


0 5 


D ±1 

^4±2,±i ^4±i,o ^4±2,o 


(S24) 


£*±2 D±i 


D± 2 


0 5 


which were also collected into Eq. (4) in the main text. The transmission coefficients, Eq.(|S23|), supplemented by 


Eq.(S24) have a physically transparent form if th ey describe propagating waves, i.e. waves with all E n > \k y \. In 

this case we can identify D~ x = tn^ [see Eq.(Sll)] and, reading the resulting expressions from right-to-left, we can 
distinguish the transmission processes depicted in Fig. 1(b) in the main text [except that E\ ^ E 5 in the process (3), 
according to our assumption]. 
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B. Close-to-resonance transmission 


Now we will focus on the resonances associated with the case when the energy of one of the n = ± 1 side-bands hits 
the bound state, E±i = E and the corresponding channel is closed. They are observed as zeros in T ±i and maxima 
in T± 2 , dispersing with the incidence angle cp (see Fig. SI). For definiteness we will consider the resonance condition 
for n = 1, but this analysis is straightforward to repeat for n = — 1. So, the resonance condition reads 


D\ (k y , E r ) — 0. 


We expand the D\ coefficient in Eq.(S21) around the resonance energy assuming 

E = Ef) — T SE = E r T SE , \SE\ {U, |*^|}, 
A K 

D\{ky , U) 


1 * 2/1 sin 2 Z 0 


(S25) 


(S26) 


Evaluating all other functions in Eqs. (S21) at E = E r , we solve the resulting system of equations keeping only terms 
of order 0(SE , Zf). The solution for to and reads 


to 


t2 = 


SED 2 — ZfAi^A 2 ,i\k y \ sin 2 Zq 


SED0D2 — Z\ (D2 Aq,i Ai^o + ^ 0 ^. 1 , 2 ^ 2 , 1 ) |* 2 /| sin 2 Zo 
_ ZfA 2i iA li0 \ky \ sin 2 Z 0 _ 

SED 0 D 2 — Z\ (112^-0,1^1,0 + Al^- 1 , 2 ^ 2 ,l) 1*2/1 sin 2 Zo 


(527) 

(528) 


We note that for \k y \ 0, |£o| 2 will be close to unity due to Klein tunneling [H9] and there is no resonance behavior. 

If w e con si der t he case when both the main channel n = 0 and the second side-band n = 2 are propagating, then 
Eqs.(S27)-(|S28|) can be rewritten as 


SE - ZfA h2 t^A 2A \k y \ sin 2 Z 0 


to = 


t2 = 


SE - Zl (4 0) ,40,1^1,0 + ^ 1>2 4 0) A 2>1 ) \k y \ sin 2 Z 0 
ZlA 2A A lfi t^\k y \sm 2 Z 0 


SE 


~ (4 0) Ao,iAi,o + ^i,24°^2,i) \ky\ sin 2 


.(0) 

0 5 

(S29) 

.(0) 

'2 5 

(S30) 


with the short hand notation t q = t^°\E r ) and t 2 = t^°\E r + 211). If we analyze these expressions we can see the 
following: 

• if we compute the corresponding transmission probabilities, Tq = |to | 2 , and X 2 = |^ 21 2 , we clearly see that Tq has 
a Fano-type resonance shape, while X 2 is of Breit-Wigner type with the width of the resonances oc Z\\k y \ sin 2 Zo . 

• exactly at the resonance, SE 0, both £ 0 and £2 have finite values independent of Zi due to constructive 
interference between the first and the second side-bands. 


In Fig. S2 we compare the approximate solution we have f ound with the exact numerical calculation [see inset of 
Fig.2(b) in the main text]. We clearly see that Eqs.(S29)-(S30) correctly describe all the essential features of the 
transmission probabilities discussed above. 


Finally, in a more strict expansion of all functions in Eq. (S21) to linear order in SE , more cumbersome expressions 


are obtained, but the above conclusions will not change, as also supported by the good agreement between the black 


and red lines in Fig. S2 


IV. SCATTERING FIELD THEORY OF AC CURRENT 


In this section we briefly describe the method we used to compute AC electric current. The theory below is valid 
as soon as a single-particle approach is justified, i.e. when particle-particle interactions can be neglected. Without 
loss of generality we assume particles incident on the barrier from the contact a [e.g. the source contact, see Fig. 1(a) 
in the main text]. Using the scattering basis, Eq.(S6), found above we construct a field operator 


% 


+00 

w/ k ’' / 


|£|>IM 


dE 

\Z2tt 


—iEt 


['7cK,in(* 7 /’ -*'0^—>-(*U * 2 /’ ■*') To;,out(* 2 /’ -^)VV- (*U *2/’ -^0] > 


(S31) 
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FIG. S2. Energy dependence of transmission probabilities To and T 2 for incidence angle ip = 7r/9, Zo = 0.47T and Z\ — 0.45. 
The blue dotted line is the result for To when neglecting scattering to the second sideband at E<i. In this case, the Fano 
resonance is fully developed (peak at unit transmission and dip at zero transmission). 


in the local coordinate system of the contact, where 7 a ,in/out(^y? E) are the corresponding annihilation operators for 
the incoming/outgoing particles, which satisfy 

{% M (k y ,E)^l. n (k' y ,E')} = 6 a , p 5(k y - k' y )5(E - E 1 ), (S32) 

{'la,in(ky, E), E )} = jybn (ky, E ), O^in (k y , E ) j- =0. 

According to the scattering theory the outgoing operator 7 a ,out {k y ,E) is, via a scattering matrix, related to the 
incoming one. For our case of an AC barrier and static contacts this relation reads 

7 c*,out {ky i E') E E S a /3(k y ', E, E n )qfp t i n (ky, E n ), (S33) 

/3 n, propag. 

where we restrict the sum over side-bands to propagating waves only, which is equivalent to setting the scattering 
matrix elements to zero if an incoming/outgoing wave is evanescent. Then we construct the current operator defined 
by the standard expression [^3] 


I a (x,t) = e 


J dy 'HC 


x,y,t)a x $ a (x,y,t), 


(S34) 


where e is the electron charge. Note that 8{k y — k' y ) in Eq.(S32) must be understood in a sense of a Kronecker symbol 
meaning that we use Born-von Karman periodic boundary conditions in the y-direction. It means that there is a 
correspondence 


5{k y -k’ y )= J fdye i W-Wv = 8 n , m , 

— oo ^ 0 

+oo 

^ / dky 

n. J 


27r 

L v 


y —oo 


(535) 

(536) 
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To obtain an observable quantity I a (x,t) we compute a statistical average of Eq.( |S34| ) with the help of 

{l ] aM {k y ,E)^ M {k y ,E')) = 5 a ,p5(E-E')f a (E), 

where f a (E) is a Fermi-Dirac distribution in the contact a. The resulting expression has the form 


+oo 


I a (x,t)= Y, e~ inm I a ,n(x), I a ,-n(x) = /*,„(*), 


(S37) 


(S38) 


where 


+ OO /- 

^ J* dk y j dFx d n Q f a (E^ 

— I ZT 1 1 \ I U I 


\E\>\ky\ 


+ == = e-^y^+^y^S^iky, E n , E)f a (E) 

2y/v(k y ,E)v(k y ,E n ) 

+ ==^ = e <["-(fc«.g-n)+"-(fc„g)]» [S aa (k y -E_ n , £)] f f a (E) 

2V v ( k y’E-n)v(k y ,E) 

, V' rfVhtlYl + J[Kx(ky,E)-K x (k y ,E n )]x 

2 1 /v(A; y ,£)v(A; y ,£ n ) 

x [S , Q/3 (fc ?/ ;-B,£’ m )] t S' Q/ 3(A:j / ;£’ n ,£; m )/ / 3(£; m ) >. 


(S39) 


Using unitarity of the scattering matrix [qlO], 

E n , E m )]^ S a j(k y ; E n , F 7 ) — ^/ 3 ,7^777,05 

a n 

YY s ^ e 777 5 En) [*^Q;^(^7/5 E : -E77,)] 7^777,0 5 

/3 77 


we can rewrite Eq.(S39) in the following form 

+00 


-hoo ✓ 

, n (x) = e[dk y [ dE\ E") e -i[ Kx (k y ,E) +Kx (k y ,E n)]xSaa{k En ^E)f a {E) 

_J J [2 ^/v(k y ,E)v(k y ,E n ) 


l^l>|fevl 


Tj ( ky 5 E— 77) + TJ (ky , E) j\ Kx (k V 'E- n ) + K x (k v ,E)]x 


2y/v{ky,E_ n )v{k y ,E) 


77 (k y ,E m ) + T}(ky, E n + m ) ^j\ Kx (k v ,E rn )-K x (k v ,E ri + rn )]3L 

p 777=-00 2y / v(/c ?/ , E m )v(k y , F'n+777) 
x [S'a/3(^;E , m ,E 1 )] t 5 aj g(fe y ;E n+m ,£?) [fp(E) - f a (E m )\ >. 


(540) 

(541) 


v>jfc)JaJ [5 aa (*y;S-n,^)] T /a(^) (S42) 


In contrast with the usua l Biittiker theory [^5], one cannot in general neglect the energy dependence of K x (k y ,E) 
and y(k y ,E m ) in Eq.(S42), because the Fermi energy Ep in graphene can be tuned to the Dirac point. On the other 
hand, if we keep the first two terms on the rhs of Eq.(S42), we see that the AC current is formally determined by the 


full Fermi sea rather than states close to the Fermi surface only. 


V. AC DIFFERENTIAL CONDUCTANCE 

In this section we present formulas that we use to compute AC conductance for different side-bands in the main text. 
We assume that our system [see Fig. 1(a) in the main text] is at low temperature and compute a linear differential 
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(a) 0 . 45 -j 


\ \ 

/ / 


-0.45 ■] f p-0-175 






(d) 0.45-1 


\ \ 

/ / 



FIG. S3. Source-drain linear conductances for side-band currents with n = ±1 and n — d=2 in the presence of ac drive on 
the top gate. Upper panels: impact angle resolved average conductances G n (EF,(p)- Lower panels: angle integrated real and 
imaginary parts of average conductances G u (Ef)- 


conductance with respect to the source-drain bias voltage Vs, 


f a (E) = f(E - eV a ), 


Gu{E f ) = 


dE 

dlD,n(% = 0+ Vs) 


9V S 


(S43) 


V s ->0 


Note that in principle the current, Eq. (S38), is a function of coordinate and we choose the point x = 0 + in our 


calculations. If we use the results of the previous section we obtain 


G r 


{Ef) = T f dk « £ y7=====4(^,^n +m (^^) 

^ J m=—oo ^ V v (ky, E m )v(ky, E nF m) 


(S44) 


E=Ep 


where we have restored h to obtain the well-known conductance unit. This for mula was used in Fig. (3) of the main 
text. Finally, in Fig. (S3) we present the results obtained with the help of Eq. (S44) for side-bands with n = ±1,±2. 
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